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DEPENDENCE  OF  REMOTE  TEMPERATURE  RETRIEVAL  ON  ATMOSPHERIC 
TRANSMITTANCE  ACCURACY 


I 


J.I.F.  King 


Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  AFB,  MA  01731-3010 


In  the  remote  temperature  sensing  problem  the  earth-atmospheric  system  is  scanned  across  the 
far-infrared  emission  bands  of  CO2  and  HjO.  The  Planck  radiation  profile  is  degraded  in  its 
upward  passage  by  scattering,  absorption,  and  re-emission  interactions  with  the  transmitting 
medium.  The  effects  of  inaccurate  atmospheric  transmittances  on  the  temperature  inferences 
from  the  satellite  data  will  be  demonstrated  and  discussed. 
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MAPPING  ATMOSPHERIC  WATER  VAPOR  AND  THE  INVERSION  OF  SPECTRAL 
RADIANCE  TO  APPARENT  REFLECTANCE  WITH  MODTRAN2  AND  DATA 
MEASURED  BY  THE  AIRBORNE  VISIBLE-INFRARED  IMAGING  SPECTROMETER 
(AVmiS) 


R.O.  Green,  J.E.  Conel,  T.G.  Chrien 


Jet  Propulsion  Laboratory 
California  Institute  of  Technology 


The  Airborne  Visible-Infrared  Imaging  Spectrometer  (AVIRIS)  is  a  NASA-sponsored  Earth¬ 
looking  imaging  spectrometer  designed,  built  and  operated  by  the  Jet  Propulsion  Laboratory. 
AVIRIS  acquires  flight  data  from  the  Q-bay  of  a  NASA  ER-2  at  20  km  altitude  that  is  operated 
from  the  Ames  Research  Center.  This  imaging  spectrometer  measures  the  total  upwelling 
spectral  radiance  from  400  to  2500  nm  in  the  electromagnetic  spectrum  through  224  channels 
at  10  nm  spectral  intervals.  Data  are  acquired  as  11  km  by  up  to  100  km  images  with  20  m  by 
20  m  spatial  resolution.  The  spectral,  radiometric  and  geometric  characteristics  of  AVIRIS  are 
calibrated  in  the  laboratory.  These  characteristics  are  validated  through  inflight  calibration 
experiments  that  are  carried  out  each  year.  MODTRAN2  is  used  for  mapping  atmospheric  water 
vapor  and  for  the  inversion  of  spectral  radiance  to  apparent  surface  reflectance  with  data 
measured  by  AVIRIS. 
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OVERVIEW 
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AVIRIS  SPECTRAL  RANGE  AND  RESOLUTION 
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Calibration:  atmosphere,  satellites,  aircraft  systems 
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Atmospheric  Transmission  and  Plant  Reflectance 


Wavelength  (nm) 


Reflectance  for  NLLS  fit  with  MODTRAN2 
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Wavelength  (nm) 


AVIRIS:  MEASUREMENT,  CALIBRATION 
AND  CALCULATED  REFLECTANCE 
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Wavelength  (nm) 


AVIRIS  Apparent  Reflectance,  Jasper  Ridge:  Forest  Target 
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AVIRIS  Apparent  Reflectance,  Jasper  Ridge:  Forest  Target  (incorrect  water  vapor) 
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Wavelength  (nm) 


AVIRIS  Apparent  Reflectance,  Jasper  Ridge:  Forest  Target  (incorrect  path  radiance) 
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THE  PRESENCE  OF  EARTH  ATMOSPHERIC 
ABSORPTION  BANDS  IN  THE  LOWTRAN7  SOLAR 
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INTRODUCTION 


The  ATMOS  Experiment  during  Spacelab  3 
Mission  in  the  Spring  of  1985. 

Reduction  of  ATMOS  spectra  above  the  earth 
atmosphere  to  obtain  a  transmittance 
spectrum  of  the  solar  atmosphere 

Comparison  of  the  ATMOS  spectrum  with 
LOWTRAN7  and  other  solar  irradiance  curves 

Problems  with  solar  irradiance  measurements 
in  late  1960's  and  early  1970’s. 

Improved  measurements  of  solar  irradiance 
spectra  in  the  visible  and  near-IR  from  satellite 
platforms  are  necessary  in  the  future. 


THE  ATMOS  EXPERIMENT 
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SPECTRAL  BACKGROUND  LEVEL 


ATMOS  SOLAR  TRANSMITTANCE  SPECTRUM 
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Improved  measurements  of  solar  irradiances  in  the  visible  and 
near-IR  region  with  accurate  radiometric  calibrations  from  satellite 
platforms  are  necessary  in  the  future. 


THE  IMPACT  OF  THIN  CIRRUS  CLOUDS  ON  TERRAIN  REMOTE  SENSING 


W.M.  Cornette,  J.G.  Shanks 


Photon  Research  Associates,  Inc. 

10350  N.  Torrey  Pines  Road,  Suite  300 
La  Jolla,  CA  92037-1020 


The  existence  of  thin,  or  sub-visual,  cirrus  clouds  within  the  field-of-view  of  a  satellite  sensor 
can  alter  the  sensor’s  perception  of  the  terrain  radiance.  Since  the  existence  of  sub-visual  cirrus 
is  frequently  unknown,  this  alteration  of  the  terrain  radiance  can  lead  to  erroneous  conclusions 
from  analysis  using  this  data.  Even  if  the  existence  of  thin  cirrus  is  known  to  exist,  the  cirrus 
parameters  (e.g.,  cloud  optical  depth,  cloud  altitude  and  thickness,  particle  size  distribution)  are 
not  fully  known.  This  paper  will  quickly  review  cirrus  cloud  parameters  and  how  these 
parameters  are  used  in  MODTRAN  to  evaluate  the  impact  of  thin  cirrus  clouds  on  remote 
sensing.  The  difference  between  actual  and  predicted  surface  temperatures  in  the  presence  of 
cirrus  will  be  presented.  A  couple  of  numerical  problems  with  MODTRAN  that  were 
encountered  during  this  study  will  be  discussed. 
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fi  042  411 


1  < -  277  -  T  1  apparent  296 . 

Fig-A.  Apparent  Temperature  Image,  Thin  Cirrus  Over  Desert 


-k  Spectral  Band:  [11.4,  12.41  (um) ,  (=  AVHRR  Band  5) 
1024  X  1024  ©  0.12  (kn)  (from  LAHDSAT),  Hist  Eq.  Scale 
^  SSGM  R  5.0  /  Modified  CLDSIM,  GBIESSIS  Scenes,  3/31/93 
>  Joe  Shanks,  Photon  Research  Associates,  Inc.  ,  La  Jolla 
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Thermal  Load 


TEMPERATURE  INVERSION  EQUATION 
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CIRRUS  CLOUD  PARAMETERS 
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Phase  Function  and  Asymmetry  Parameter 


AHENUATION  COEFFICIENT 
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cinnus  THICKNESS  l  (Ki..)  temperature  (K) 
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APPLICATION  OF  LOWTRAN-7  TO  AVHRR  THERMAL  DATA  IN  FIFE 


T.  Schmugge,  P.  Bougarel  M.  Sugita 


W.  Brutsaert 


USDA  Hydrology  Lab  Env.  Res.  Ctr.  Cornell  University 

Beltsville,  MD  Univ.  of  Tsukuba 

Japan 


LOWTRAN-7  was  used  with  near  coincident  radio  soundings  of  the  atmosphere  to  estimate 
surface  brightness  temperatures  from  the  NOAA-9  AVHRR  data.  The  data  were  obtained  for 
8  days  during  the  First  ISLSCP  Field  Experiment  (FIFE)  conducted  in  central  Kansas  between 
the  end  of  June  and  October  of  1987.  In  general  the  channel  4  results  were  up  to  2”  K  warmer 
than  those  from  channel  5  indicating  that  we  were  undercorrecting  for  the  atmospheric  water 
vapor.  However  on  one  day  the  reverse  was  true  and  the  channel  5  result  was  higher.  In  this 
case  we  were  able  to  get  agreement  by  reducing  water  vapor  content  to  80%  of  its  measpured 
value.  The  results  were  compared  also  with  various  split  window  approaches  for  estimating 
surface  temperature  and  with  ground  based  broad  band  radiometers  with  AVHRR  results  being 
generally  warmer. 
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Application  of  LOWTRAN-7  to  AVHRR  Thermal  Data  In  FIFE 
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ABSTRACT  for  Annual  Review  Conference  on  Atmospheric  Transmission 
Models 


Application  of  LOWTRAN-7  to  AVHRR  Thermal  Data  In  FIFE 
T.  Schmugge  and  P.  Bougarel,  USDA  Hydrology  Lab,  Beltsville  MD 
M.  Sugita,  Environmental  Research  Center,  Univ.  of  Tsukuba,  Japan 
W.  Brutsaert,  Cornell  University 

LOWTRAN-7  was  used  with  near  coincident  radio  soundings  of  the 
atmosphere  to  estimate  surface  brightness  temperatures  from  the 
NOAA-9  avhrr  data.  The  data  were  obtained  for  8  days  during  the 
First  ISLSCP  Field  Experiment  (FIFE)  conducted  in  central  Kansas 
between  the  end  of  June  and  October  of  1987.  In  general  the 
channel  4  results  were  up  to  2  K  warmer  than  those  from  channel  5 
indicating  that  we  were  undercorrecting  for  the  atmospheric  water 
vapor.  However  on  one  day  the  reverse  was  true  and  the  channel  5 
result  was  higher.  In  this  case  we  were  able  to  get  agreement  by 
reducing  water  vapor  content  to  80%  of  its  measured  value.  The 
results  were  compared  also  with  various  split  window  approaches 
for  estimating  surface  temperature  and  with  ground  based  broad 
band  radiometers  with  avhrr  results  being  generally  warmer. 
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®  SP-Super  PAM  j 
®  SD-Super  DCpJ 


•  B-Bowen  Ratio  Rux  Measurement 
X  E-Eddy  Correlation  Rux  Measurement 
— —  Site  Boundary 

iiiiiii  1-70 

===  R-177 


Figure  2.2b;  Map  of  FIFE  sites  showing  flux  stations  and  AMS  as  actually  located 
during  FIF^  987.  Stations  are  marked  by  station  number:  even 
numbers  for  flux  stations,  odd  numbers  for  AMS. 

Note  (i)  No  DCFs  displayed,  as  data  is  currently  inaccessible 

(ii)  Stations  40  (flux)  and  31  (AMS)  moved  from  positions  as  marked 
in  EXPLAN 

(iii)  25  (east)  and  29  (north)  operating  in  IFC-4  only. 

See  also  Table  2.2b  and  Figure/Table  2.4b  in  EXPLAN 

2-  5 


73 


r-^  '  '  '  '  'CO 

-o  o  o  o  o  o 

o  03  CD  CM 


a3iillMSNViJi  % 


74 


OATS 


x: 

o 

CD 

O 

o. 

Q. 

CD 


O 

U 

C 


Q. 

C/3 

DC 

DC 

UJ 


C/D 

c 

0 


0 

V- 

D 

0 

0 

0 

E 

■0 

c 

D 

O 

k- 

OD 


C 

o 

‘v» 

0 

Q. 

E 

o 

o 

> 

0 

Q. 

0 

> 

0 

Q 


Lf) 

# 

CD 

csj 


.  =  H 
0  ® 
0  CO 

CJ 

0 

D  + 

"D 

0 

+=>  ' 

0  II 

0 

D)  H 
0  =3 

>  ^co 

O 

UL 


in 


CM 


h" 


OQ 

+ 


’o 

CO 

0 

II 

0 

'■4-” 

W. 

0 

D 

C/5 

k- 

0 

UL 

h“ 

0 

0 

■M 

’0 

"O 

c 

D 

O 

u. 

DD 

0 

JZ 


O 

0 

0 


0 

0 

O 

0 

H” 

l_ 

o 

0 

"^3 

0 

.1^ 

0 

"M 

0 

DD 

0 

> 

"O 

0 


D 

0 

0 

0 

0 


CM 

1 


ch* 


1 


s 

o 

L 

::i 

s: 

di 


X 

Od 


j. 


76 


315  320 


82 


84 


c 

0 


D 

O 

c 

o 

o 


0 


>; 

JD 

0 

C 

o 

w 

0 

0 

I  “ 

O  0 

o  w 

4-^  O 

52  E 

^  0 
0 

Q.  0 
ClSZ 
0 


0 

c 

o 

■ 

+-» 

0 

> 

!_ 

0 

0 

o 

■D 

c 

D 

O 

05 


o 

CO 


0 

■M 

+-> 

0 

0  -Q 

-52  > 

05  +-» 

c  x: 

0  05 


^  E^  0 


c  o 

0 

0) 

0 

J_  0 

0 

o 

z 

c 

o 

o  .E 

0 

C/) 

^  c 

0 

H- 

D 

c 

■ 

_J 

0  o 

Q 

o 

»-  0 

z 

1  ° 

' 

o 

u 

0 

c 

0 

> 

0 


c 

0 

k- 

4-> 

T3 

O 


85 


AN  URGENT  NEED  OF  VALIDATING  WATER  VAPOR  ABSORPTION  COEFFICIENTS 
FOR  THE  DEVELOPMENT  OF  EOS’S  EARTH  SURFACE  TEMPERATURE 
ALGORITHMS 


Z.  Wan,  J.  Dozier 


Computer  Systems  Laboratory 
Center  for  Remote  Sensing  &  Env.  Optics 
Univ.  of  California 
Santa  Barbara,  CA  93106 


In  order  to  remotely  measure  surface  temperature  from  space  to  an  accuracy  of  O.^K  for 
oceans  and  1“K  for  land,  as  required  by  NASA’s  Earth  Observing  System  (EOS),  radiative 
transfer  simulations  should  be  accurate  to  tenths  of  a  percent.  Through  "exponential-sum-fit" 
tables,  the  atmospheric  transmission  models  in  LOWTRAN  and  MODTRAN  have  been 
incorporated  into  accurate  radiative  transfer  models  based  on  adding/doubling  or  discrete- 
ordinate  methods.  Simulated  cross-comparisons  show  that  approximations  used  in  LOWTRAN 
and  MODTRAN  can  cause  AVHRR  band  radiances  to  differ  by  as  much  as  1%,  and  the 
uncertainties  caused  by  absorption  by  the  water  vapor  continuum  are  larger.  In  combination 
with  the  watepr  vapor  band  absorption  coefficients  in  MODTRAN,  with  a  finer  spectral 
resolution  and  dependence  on  temperature  and  pressure,  an  exponential  form  of  water  vapor 
continuum  absorption  with  a  factor  of  1.157  on  the  coefficient  at  296“K  used  in  LOWTRAN7 
gave  a  better  agreement  between  AVHRR  data  and  simulations  that  were  based  on  radiometric 
SST  data  and  measured  atmospheric  temperature  and  humidity  profiles.  Accurate  measurements 
and  validations  of  water  vapor  absorption  in  a  wide  temperature  range  at  a  moderate  spectral 
resolution,  5  cm'*,  are  recommended. 
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need  to  check  various  approximations  used  in  radiative  transfer 
modeis, 

need  more  accurate  water  vapor  absorption  coefficients. 


TABLE  J. 
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ATRAD 

—  deals  with  accurate  multiple  scattering  based  on  the  interaction 
principle  and  adding/doubling  method 
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•  Comparison  between  MODTRAN  and  LOWTRAN 

TABLE  2.  The  band  transmission  and  brightness  temperature  of  MODIS 

and  ASTER  thermal  bands  for  a  land  surface  with  e=0.98  and  at 


0 

O 

0 


13 

C/3 


£  03 

0  E 

JZ  ■> 

s-l 

o  > 


cd  ^ 
o 


CO 

h- 

I 

o 

5 


o 


o 


N. 

h" 

I 

o 

S 
I — 


C  I 
CO  g 


"d 

0 

C 

0 

E 

1 

CO 

c 

C0 

0 

"D 

LO 

< 

C 

LO 

0 

13 

0 

CD 

C 

CO 

0 

0 

"D 

fx. 

E 

C 

CO 

03 

V 

rx 

03 

(M 

CO 

11 

CM 

0 

CO  . 

.i- 

c 

CD 

h? 

■*— » 

CO 

II 

Id 

.  ^ 

'0 

°> 

o 

c 


Q 

O 


<X)^CJ>CNCON. 
^  ^  CM  CO  CM  ^ 

d  CD  CD  d  ^ 


N.  CD  CO  CM  T“  ID 
CM  to  CO  CM  CD 

d  d  d  d  d  d 


vP  vP  xP  v2  vO  vO 
^  ^  ^  ^  ^  ^ 
T—  CO  O  ^  CD 


Cd  T-  <M  T 

I  I  I  -J.  4.  ^ 


xP  vp  vp  o  vP  vP 

^  ^ 

qCDN;;^?-!^ 

CD  'f”  CO  o  d  CM 
+  4-  +  »  *  + 


O^lDOCDI^ 
^  CD  O  CM  CM 

00  CD  o 

CO 

t  I  I  I  I  I 

O^UDOCON. 

CDCOCMOI^rx 

q  q  q  ^  Q  ^ 

CO  CO  ^  00  ^  ^ 


O  CM  CO  CD  'r-  CM 
CM  CM  CM  CM  CO  CO 


^  00  O  O  CD 
N-  q  CD  00  o 

O  O  T“  *r“ 


LOcocDinr^ 

CD  CO  ^  O 

d  d  d  d  d 


CD 


Vp  xp  vP  vP 


00  00  CD  00 

d  CD  CD  CD 
+  +  -f  + 


xp  xp  "xp  xP 

rx.  ^  q  q  cb 

CM  d  V  T  o 


LT)  m  10  Lo  in 

CM  CD  CD 

^  q  q  d  ^ 

00  00  CD  'T-  ^ 

I  I  1  I  I 

Ln  in  m  m  m 

c\j  x  C\i  C\j  O) 

"Z  ^  o  o 

cn  CO  C30  -r-  X- 


cr 

iflo-^CVJCOTf 

^  -r-  X-  X-  ^  -r- 

< 


94 


•  Factor  Analysis 

TABLE  3.  Comparison  of  molecular  absorptions  in  different  models  with 
MODTRAN. 
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Note:  (+,-)  means  that  it  is  larger  in  NOAA-7  AVHRR  channel  4  but  smaller  in  channel  5. 
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TABLE  4.  Quantitative  comparisons  between  different  models  in  NOAA7 
AVHRR  band  brightness  temperature  and  MCSST  values 
simulated  for  tropical  atmosphere,  surface  visibility  23  km  at 
0.55  jim,  8=0.98,  Tss=Ta/r=  299.7  °K,  viewing  at  11.4°. 
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__  (°K) 

Ts 
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74-7. 

(°K)  " 
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(°K)  -  SST 
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MODTRAN 

295.434 

293.509 

1.925 

300.57  +0.87 
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ATRAD-MOD1 

295.175 
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1.415 

298.99  -0.71 

2’ 

ATRAD-MOD2 
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299.37  -0.33 
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ATRAD-LOW1 

294.674 

293.207 

1.467 

298.60  -1.10 
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H2O  band  difference 
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no  k-distribution 
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Comparison  with  radiative  transfer  simulations 


TABLE  6,  RT  simulation  results  compared  with  SST  measurements. 
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The  sea-surface  emissivity  model  (Masuda  et  al.,  1988)  has  been  used. 

"exptbl{lowtran7)"  means  the  exponential-sum-fitting  table  from  the  transmission  values  used  in 
L0WTRAN7. 

"exp.  form"  means  the  exponential  form  of  the  H2O  self-broading  cofficient  In  Cg  =  Q/T+constant 
based  on  the  coefficient  at  296  “K  in  LOWTRAN7,  and  data  of  Burch  and  Gryvnak  (1980). 


STAUS  OF  VALIDATION  OF  WATER  VAPOR  ABSORPTION  COEFFICIENTS 
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Figure  4.  Comparison  between  atmospheric  transmission  dam  and  RT  simulations. 
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ANALYSIS  OF  CLOUD-TOP  HEIGHT  AND  RELATED  CLOUD  PARAMETERS 
FROM  SATELLITES  USING  THE  Oj  A  AND  B  BANDS 


A.  Kuze,  K.  Chance 


Harvard-Smithsonian  Ctr. 
for  Astrophysics 

Atomic  &  Molecular  Physics  Div. 
60  Garden  Street 
Cambridge,  MA  02138 


Cloud  height  and  cloud  coverage  detection  are  important  for  total  ozone  retrieval  using 
ultraviolet  scattered  light.  Use  of  the  O2  A  and  B  bands,  around  762  and  688  nm,  make  it 
possible  to  detect  both  cloud  top  height  and  cloud  percentage.  The  measured  values  of  a  space 
borne  high  resolution  spectrometer  are  convolutions  of  the  instrument  slit  function  and  the 
atmospheric  transmittance  between  cloud  top  and  satellite.  Optical  thicknesses  between  the 
satellite  orbit  and  each  height  are  calculated  with  high  accuracy  using  FASCODE3P.  Cloud 
parameters  are  determined  by  least-squares  fitting.  A  grid-search  method  is  used  to  search  the 
parameter  space  of  cloud  top  height  and  percentage  to  minimize  the  variance.  For  this  search, 
nonlinearity  of  atmospheric  transmittance  is  important.  Using  the  above-mentioned  method, 
operational  cloud  detection  is  possible  with  minimal  computation  time.  Measurement  of  clouds 
and  atmospheric  trace  gases  in  the  same  IFOV  are  also  possible. 
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Analysis  of 

Cloud-top  Height  and  Related  Cloud  Parameters 
from  Satellites  Using  the  O2  A  and  B  Bands 


June  8,  1993 


AKIHIKO  KUZE  AND  KELLY  V.  CHANCE 


Harvard-Smithsonian  Center  for  Astrophysics 


IFOV  of  Satellite  Borne  High  IFOV  of  Radiometer 

Resolution  Spectrometer 
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Ozone  Profile  and  Cloud  Effect  on  Measurement 
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GOME  and  SCIAMACHY  Project 

Satellite  Borne  High  Resolution  Spectrometer 
Orbit  :  Sun-Synchronous 
Altitude  :  800  km 

Nadir  Looking 

Solar  Backscattering  Radiance  Observation 
IFOV  :  320  km  x  40  km 
Spectral  Coverage  :  240  -  790  nm 
Spectral  Resolution  :  0,1  -  0.2  nm 
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Radiative  Transfer  in  the  O2  Absorption  Re 
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aij  :  the  cloud  top  reflectivity 
Pj  :the  earth’s  surface  diffusive  albedo 
r(zy,  h)  i  the  optical  depth  between  cloud  top  and  satellite 
(proportional  to  total  O2  column  amount  al 
fji^)  •  ^he  slit  function  of  channel 
F(i')  :  the  solar  spectrum 
ri  :  the  coverage  of  type  i  cloud 
F  :  the  too  heie:ht  of  tvoe  i  cloud 


Ai/  :  FWHM  of  the  instrument 
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02  A  and  B  band 
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Algorithm 


1 


Q-value  table 
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Height  and  Coverage  Detection 
Optical  Depth  Effect 
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no 


Height  and  Coverage  Detection 
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Low  Height  &  Low  Coverage  Low  Height  &  H  igh  Coverag 


(Uoiid  top  height  (kin) 


Accuracy  Requirements 
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Error  Analysis 

Measurement  Error 

Model  Calculation  Error 
Scattering 
Earth  Albedo  Data 
Spectral  Band  Database 

Intensity  and  Line  Width 
Slit  Function 
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Conclusion 


Characteristic  of  this  Method 

Both  Cloud  Top  Height  and  Coverage  Detection 
with  Minimal  Computation  Time 
in  Moderate  Accuracy 


Further  Studies 

Spectral  Band  Database 

Earth  Surface  Albedo  Information 
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THE  USE  OF  SPACE  BASED  REMOTE  SENSING  FOR  ESTIMATION  OF  THE 
METHANE  MIXING  RATIO  IN  THE  MIXING  LAYER 


P.  Ashcroft 


Dept,  of  Engineering  &  Public  Policy 
Carnegie  Mellon  University 
Pittsburg,  PA  15213 


This  investigation  explores  the  limits  of  space  based  remote  sensing  for  methane 
characterization,  and  the  extent  to  which  those  limits  are  affected  by  uncertainty  in  other 
atmospheric  attributes.  Central  to  this  analysis  is  the  fact  that  the  attribute  considered, 
(enhancement  of  the  methane  mixing  ratio  in  the  mixing  layer),  is  not  directly  observable,  but 
must  be  obtained  through  inversion  of  the  spatially  integrated  signal  received  by  the  instrument. 

MODTRAN  is  used  to  simulate  atmospheric  transmission  in  the  neighborhood  of  the  3.3^ 
methane  absorption  band.  This  signal  is  then  inverted  to  determine  the  importance  of  ancillary 
information,  (e.g.,  the  vertical  profile  of  HjO  and  temperature),  for  methane  retrieval.  Thus, 
this  research  provides  quantitative  description  of  the  importance  of  simultaneous  observation  of 
methane  with  other  atmospheric  conditions,  and  the  relative  utility  of  space  based  observation 
relative  to  other  measurement  methods. 
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Remote  Sensing  Policy 
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Examples  of  Remote  Sensing 


Case  Study:  Atmospheric 
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Specific  Question 
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analytic  Approach 
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•  This  is  the  "Best  Case"  for  methane  retrieval  subject 
to  uncertainty  in  the  water  column,  but  assuming  that  there 
are  no  other  sources  of  uncertainty. 


assumptions  for  this  Case  Study 


Linearized  Inversion  of  Total 
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FJ^) 


Total  Transmittance 
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Wavenumber  (cm-1) 


Gross  Water  and  Methane 
Signal  Comparison 
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Weighting  Parameter 
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Error  in  Methane  Retrievals  as  a  Result  of  1%  Error  in  Water  Profile 
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Weigliling  Parameter 
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•  Separating  the  water  measuring  instruments  from  the 
methane  measuring  instruments  may  critically  threaten 
the  ability  to  measure  methane  at  all. 


Future  work 
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Open  Questions 
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•  What  problems  should  I  anticipate  in  Implementation? 


A  CORRELATED  K-DISTRroUTION  MODEL  OF  THE  ATMOSPHERIC  HEATING 
RATES  FOR  OVERLAPPING  SPECTRA  OF  COj/H^O  AND  CH^/NjO 

% 


ALLEN  S.  GROSSMAN 
KEITH  E.  GRANT 


Global  Climate  Research  Division 
Lawrence  Livermore  National  Laboratory 


The  problem  of  overlapping  the  spectra  of  COj  and  HjO  in  the  0-2500  cm’*  wavenumber 
region  has  been  considered  using  a  correlated  k-distiibution  model  for  the  transmission  in  the 
atmosphere  between  0  and  60  km.  Individual  k-distributions  for  the  two  gases  have  been 
combined  using  an  overlapping  algorithm  to  produce  k-distributions  for  the  mixture  of  COj  plus 
HjO.  A  prototype  radiative  transfer  model  has  been  used  to  calculate  atmospheric  fluxes  and 
heating  rates.  The  heating  rates  for  the  overlapped  mixture  agree  to  better  than  eight  percent  with 
those  calculated  for  a  single  pre-combined  mixture  of  CO2  and  HjO. 
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TRANSMISSION  MODEL 
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The  transmission  expression 
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lonotonic  nature  of  the  k-distriDutions  can  provide  sunicient 
utational  economy  in  transmission  calculations  to  warrant  such 
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heating  rate(deg  C/day) 
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H^O  and  CO2  model  parameters 
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0-2500  cm  ,  25  cm  intervals 
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ALTITUDE(KM) 
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OVERLAP  MODEL 
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The  derivation  of  Eq.  4  can  be  accomplished  as  follows.  The  form  of  Eq.  4  in  terms  of  the 
distribution  functions  F2(  T2),  and  the  general  expression,  Eq.  5,  is 

G,2(^)  =  j7'''"l(’"-'2)f2(’-2)<iM5"- 

From  Eq.  5  the  expression  dG2  ( ^2)  is 

dG2{T2)  =  F2{'t2)dt2. 


and 

Gi2(T)  =  0V(l^,'r2)^^2^^=£ 

Using  the  Iterated  Integral  Theorum  for  interchanging  the  integration  Umits,  the  following 
expression  can  be  obtained; 

72^72^/7' =£  ,t2)dt'  dT2 

F,{t'-t2)ddF2{T:2)dT2.  (S) 

Jo  Jr2 

Examination  of  the  term 

Fi{d-Z2)dt', 

in  Eq.  8,  by  itself,  reveals  that  the  transformation 

u  =  't'—T2, 
du  =dt' 

u  =  r-T2’  d=t 
U  =  0,  t'='V2 


can  be  made  ( T2  is  treated  as  a  constant).  Thus  Eq.  9  becomes 

'  Fi{u)du  =  Gi{u)  =  Gi{t-'V2) 

JO 

and 

j02('^2M^2  G^{'!:-t2)  =  \^Gi{T-T:2)dG{t2) 

=  Gi2(t). 


>- 
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Graphical  representation  of  the  possible  cases  for  integrating  dG12  =  G1(tau-tau2)dG2(tau2) 
over  the  interval  0  to  tau.  This  yields  a  single  pair  of  values  (G1 2(tau),  tau)  for  overlapping 
two  gases  in  the  correlated-k  formulation  of  gaseous  absorption  bands.  The  lines  of  constant 
tau  are  shown  in  tau-space.  The  actual  integration  must  be  done  in  G-space. 
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FASCODE/MODTRAN:  VALIDATION  AND  APPLICATIONS 


^G.P.  Anderson,  J.H.  Chetwynd,  F.X.  Kneizys* 

^A.  Berk,  L.S.  Bernstein,  D.C.  Robertson,  P.  Acharya 
Theriault,  '‘L.W.  Abreu,  ®S.A.  Clough,  J.L.  Moncet 


^PL/Geophysics  Directorate 
29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 


^Spectral  Sciences,  Inc. 
99  South  Bedford  St. 
Burlington,  MA  01830 


^DREV-Defence  Research  ^ONTAR  Corp. 

Establishment  Valcartier  129  University  Road 

P.O.  Box  8800  Brookline,  MA  02146 

Quebec,  Canada,  GOA  IRO 

^Atmospheric  And  Environmental  Research,  Inc. 
Cambridge,  MA  02139 


(*  Retired) 


During  the  previous  year  the  merit  of  MODTRAN2  has  been  demonstrated,  particularly  through 
comparisons  with  both  measured  interferometric  data  and  line-by-line  (LBL)  flux  divergence 
(cooling  rate)  calculations.  MODTRAN2  differs  from  MODTRAN  in  three  major  areas:  (1)  its 
spectral  data  bases  (based  on  a  two-parameter  band  model)  are  derived  from  HITRAN92  rather 
than  the  ’86  version;  (2)  the  line-of-sight  geometry  routines  have  been  modified  to  eliminate  sets 
of  errors;  and  (3)  the  radiance  algorithm  has  been  upgraded  to  respond  appropriately  to  optically 
thick  layers.  This  latter  correction  is  similar  to  those  suggested  by  Clough  et  al.  (1992)  and 
Comette  (1992)  for  FASCODE  and  LOWTRAN,  respectively. 

The  current  and  future  directions  of  the  DoD  radiance-transminance  codes  remain  robust. 
For  instance,  DOE  is  collaborating  on  a  new  effort  to  bring  about  convergence  between  the  DOE 
FASCODE-based  code  (LBLTRN)  and  FASCOD3P,  fostering  commonality  in  vectorization, 
continua,  non-LTE,  laser,  line-of-sight  and  flux  divergence  applications.  In  addition,  the  cross 
section  capabilities  (IR  CFC’s)  from  FASCODE  are  being  incorporated  into  MODTRAN2,  and 
the  development  of  joint  FASCODE/MODTRAN  inversion  algorithms  are  being  explored. 
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ATMOSPHERIC  TRANSMISSION  MODELING 
ANNUAL  REVIEW  CONFERENCE 

8  JUNE  1993 

GEOPHYSICS  DIRECTORATE/PHILLIPS  LABORATORY 


FASCODE/MODTRAN:  Validation  and  Applications 


G.P.  Anderson,  J.H.  Chetwynd,  F.X.  Kneizys 
Geophysics  Directorate/PL 

A.  Berk,  L.S.  Bernstein,  D.C.  Robertson,  P.  Acharya 

Spectral  Sciences  Inc. 

J.-M.  Theriault 

DREV/Defence  Research  Establishment  Valcartier 

L.W.  Abreu 
ONTAR,  Inc. 

S.A.  Clough,  J.-L.  Moncet 
Atmospheric  and  Environmental  Research,  Inc. 
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HISTORY 


DoD  Plan  for  Atmospheric 
Transmission  Research  and  Development 

AIR  FORCE 

o  Maintain  DoD  Standard  Atmospheric  Optical/IR  models: 

(LOWTRAN),  MODTRAN,  FASCODE,  HITRAN  Database 
o  Publish  and  Brief  Model  Updates 
o  Conduct  Annual  Tri-Service  Review 
0  Measure  and  Model  Propagation  Effects  of  the 
Free  Atmsophere 

ARMY 

o  Study  Battlefield  Conditions 

o  Develop  Models  of  Dust,  Smoke,  Chemicals,  Propagation, 
and  Diffusion  Effects 

NAVY 

0  Develop  Models  for  Marine  Environment 
o  Measure/Model  Atmospheric  Propagation 
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OPTICAL  ENVIRONMENT  DIVISION 
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PiyGP  OPTICAL  BACKGROUND  CODES 
COVERAGE  CAPABILITY 
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WAVELENGTH  (MICRONS) 


The  Problem:  RADIATIVE  TRANSFER 


1 .  The  Atmosphere  as  a  contaminant  for  E/0  Systems 

2.  The  Atmosphere  as  a  signature  source  for  natural  variability 

Solutions:  DEFINITIONS 

1 .  State  variables  (T,  p,  Cld,  Aer)  along  line-of-sight 

2.  Spectral  Characteristics  of  the  Path  Variables 

3.  Viewing  Geometry 

4.  E/0  System  Characteristics  (Spectral  Range  &  Resolution, 

Platform,  Objective) 

Solutions:  OPTIMIZATION 

1 .  Efficient  Mathematical  Algorithms  (Line-by-Line) 

2.  Accurate  Band  Model  Options 

3.  User  Friendly 

4.  Validation/Documentation 

Solutions:  DATA  ANALYSIS 

1 .  Information  Theory 

2.  Inversion  Algorithm  Development 

3.  Ground  Truth 

4.  Validation  and  Error  Estimation 
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DEFINITIONS 


Ki  =  absorption  cross  section,  related  to  molecular  properties, 
pressure  (p),  temperature  (6) 

r]i  =  column  amount  of  absorbing  (i’th)  species  =  {  n^ds 
ds  =  path  increment;  n,,  =  volumn  density 
Tj  =  optical  depth  =  k-,  rj-, 

T;  =  transmittance  =  exp  (-r) 

Tmoi  =  molecular  transmittance  =17;  =  T,  T2  T3  T4  .... 
Tj  total  transmittance  TjjjQ|-Tj.Qjj(jjjy^-T5j.gj-T^gj.QgQ] 

B(0)  =  Planck  Function  for  temperature  6 
=  Solar  (Lunar)  Source  Function 
^  =  Non-LTE  Source  Term 
W  =  Weighting  Function  =  (dTj/ds) 

OPERATIVE  EQUATION: 

=  Thermal  Radiance  =  |  B(^)  dT^  —  J  B(0)  (dTj/ds)  ds 
AND  FOR  A  SINGLE  LAYER: 

9?  =  J  B(^)  dTj  =  Bid)  [1-Tt] 
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PAGE  2 


and,  finally,  combining  thermal,  solar,  and  non-LTE  sources 
with  multiple  scattering, 

one  can  replace  the  Planck  source  function  with  a  more  general 
source  function: 

=  general  source  function  dependent  on  optical  depth  (r)  and 
viewing  geometry  =  zenith  &  azimuth  cosines) 

=  KMtt)  p  [SOLAR] 

+  [l-ojJ  B((9)  [THERMAL] 

+  ^  [N-LTE] 

+  J^s  [MUL.SCAT.] 

where: 

cOo  =  single  scattering  albedo  = 

T(^)  =  transmittance  from  top  of  atmosphere  to  layer 
p  =  scattering  phase  function 

=  m.s.  for  both  solar  and  thermal  terms  (complicated) 

and  other  terms  are  as  previously  defined!! 
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ATMOSPHERIC  STATE  VARIABLES 


MOLECULAR,  PARTICULATE,  THERMAL  PROFILES 
Some  thoughts  on  the  subject: 

5R  =  Thermal  Radiance  =  J  B(0)  dij 


1 .  First  Order  Impact  on  Radiance  Accuracy: 

6(z)  =  temperature  profile 
B(0)  =  "driving"  Planck  Function 

2.  Second  Order  Impact  on  Radiance  Accuracy: 

fi(z)  =  mixing  ratio  profiles 

d{z)  =  atmospheric  density  profile 

k(z)  =  molecular  absorption  coef  ($  and  p) 

T(z)  ==  exp(-Ilr)  =  exp(-I!/c  fjL(z)d(z)  Az) 


Given  that: 

o  the  molecular  spectroscopy  is  "in  hand",  although  always  requiring 
"investment"  in  order  to  stay  state-of-the-art, 

And: 

o  the  LARGEST  error  sources  in  the  forward  calculations  are  the 
descriptions  of  the  "state  variables",  e.g.  temperature,  pressure,  mixing 
ratio,  and  particulate  (aerosol  and  hydrometeor)  profiles. 

Then: 

o  why  work  so  hard  on  code  accuracy? 

o  why  continue  the  "investment"  in  "state-of-the-art"  spectroscopy? 

o  what  is  the  end  benefit? 
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ANSWERS  (from  a  neutral  observer!) 


*  VALIDATION  for  Signatures 

1.  targets /pollutants  against  backgrounds- 

estimate  of  variability 

estimate  of  go /no go 

estimate  of  true/false;  reduce  false  ID 

incorporate  "target  signatures " 

2.  atmospheric  specification 

local  and/or  global  "weather" 
atmospheric  contaminants  become  signatures 

3.  higher  resolution  spectroscopy  with  high  spectral 

accuracy  may  minimize  conjusion,  increasing 
signal/noise  for  both  contaminants  and  state 
specification. 

*  REAL-TIME  ANALYSIS  with  speed 

1.  Line-by-line  (LBL)  codes  and  their  derivatives  work  - 

"exact" physics  must  remain  state-of-the-art; 

Because  physicists  drive  technology  to  higher  resolution: 
lidar/laser  applications 
interferometers  at  <.OIcm-I  resolution 
remote  sensing  issues:  line  coupling,  mm-wave 
new/old/unsettled  issues: 

CFC  x-sections,  H20  continua,  C02  x  factor 
NLTE  issues 

However: 

"LBL"  is  slow,  slow,  slow!!! 

2.  Band  Models  (pragmatic,  expedient  parameterizations)  trail 
"LBL"  in  "state-of-the-art" , 

but  they  are  the  codes  of  the  NOW  and  FUTURE  for 
issues  and  analyses  with  compatible  spectral  resolution! 
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Inversion  Algorithm,  Data  Analysis 


MODTRAN  ASSUMPTIONS: 


No  overlap  between  species: 

e.g.:  r(total)  =  r(H20)  •  r(C02)  •  7(03)  •  r(CH4)  •  ... 

Line  strength  temperature  dependence  can  be  interpolated: 

e.g.:  S-Xd)  =  ^(<?„)  for  n  =  200, 225 ,250, 275, 300K 

Voigt  line  shape  is  appropriate,  based  on  "path  weighted"  p 

e.g.:  ^(p,v)  =  line  shape  function  where  oi^(p)  =  constant 

Line  contributions  falling  outside  of  bin  (Av=\  cm‘^)  but 
within  ±25  cm’^  can  be  adequately  described: 

e.g.:  "special  MODTRAN  continua" 

Two  band  model  parameters  are  sufficient: 

e.g.:  S/d  =  average  line  strength  in  bin  =  (\/Av)  USj 


n  =  ave.  #  of  lines  in  bin 
=  l/d  =  (I/Aj-)  (ly'Si)^  /  E(Si^) 


M0DTRAN2  Band  Model  Equations: 

r(total)  =  Total  Transmittance 

=  Hi  7;.  =  r(H20)  •  r(C02)  •  T(03)  •  r(CH4)  • 

and  7;.  =  {(2/Ap)  0  j  exp[-S,  u,  iE(p,v)]  dv}^ 

where: 

Av  =  I  cm^  bin 

=  (S/d)  /  (1/d)  =  absorption  coef. /density  param. 

Uj  =  path  column  am’t  of  ith  species 

X{p,v)  =  line  shape  function  where  p  =  path  pressure 

n  =  (1/d)  Ap  ~  ave.  #  of  lines  in  bin 
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CONTRAST:  Line-by-Line  vs.  Band  Model 


FASCODE  Transmittance:  (LBL) 

n(£,i)  =  =  Eoptical  depth(i)  in  layer  i 

T{v)  =  exp[-E  n(f,i)] 

or 

T{v)  =  =  T,-T,-TsTi- . 

and 

TlAv)  =  }  (Tlv)  dv  )/  Av 

MODTRAN  Transmittance:  (BM) 

Ti(Av)  =  {{HAv)  0  i  exp[-S,  Uj  ^ip,v)'\  d^}" 

T{Av)  =  n,  Ti(Av)  =  r(H20)-7IC02)-7i:03)-7(CH4)-... 
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M0DTRAN2 


New  Features: 


1 .  HITRAN92-Based  Band  Model 

2.  New  Geometry  Package  with  Consistent  Range  Determinations 

3.  New  Radiance  Algorithm  ("Linear  in  Optical  Depth") 

4.  Spectrally-Dependent  Surface  Emissivity /Albedo 

5.  Flux  Divergence  Capability 
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FREQ 

ALT 

T 

RANGE 

MEAN  TRANS 

IDIF 

1  VAR 

(cm'O 

(km) 

(K) 

(km) 

MOD(86)  FAS(86) 

(rms") 

0- 

45 

200 

500000 

.61717 

.71211 

.09494 

4.93700E-3 

500 

30 

200 

5000 

.68104 

.66882 

.01228 

5.06760E-4 

30 

300 

5000 

.71678 

.71211 

.00467 

2.698 19E-5 

15 

275 

10 

.87833 

.87800 

.00033 

2.61406E-5 

15 

225 

10 

.86676 

.86513 

.00163 

4.90820E-5 

0 

300 

.005 

.44175 

.44270 

.00095 

2.78848E-4 

0 

250 

.005 

.42247 

.41378 

.00869 

3.16083E-4 

500- 

45 

300 

500000 

.59798 

.53966 

.05832 

3.29788E-3 

1000 

30 

260 

5000 

.65704 

.61798 

.03905 

6.05537E-3 

15 

260 

500 

.57125 

.56280 

.00844 

4.04609E-3 

15 

260 

10 

.81105 

.80050 

.01055 

2.67364E-4 

0 

300 

1 

.51721 

.51560 

.00161 

4.55383E-4 

0 

260 

1 

.49087 

.48749 

.00339 

3.36086E-4 

1000- 

45 

300 

500000 

.68805 

.62673 

.06132 

1.93994E-3 

1500 

30 

260 

10000 

.63912 

.59061 

.04851 

2.39270E-3 

15 

260 

50 

.80956 

.79412 

.01544 

1.64402E-3 

0 

300 

1 

.48765 

.48703 

.00062 

2.78078E-5 

0 

260 

.1 

.68762 

.68686 

.00076 

3.01251E-5 

1500- 

45 

260 

500000 

.75053 

.72518 

.02534 

4.86761E-4 

2000 

30 

300 

10000 

.75174 

.73917 

.01257 

2.25932E-4 

30 

260 

5000 

.81098 

.79915 

.01183 

2.04952E-4 

15 

260 

500 

.55166 

.55034 

.00132 

6.32165E-5 

0 

300 

.0 

.33631 

.33473 

.00158 

1.10876E-4 

2000- 

45 

260 

5000 

.78488 

.74690 

.03798 

1.68926E-3 

2500 

30 

300 

100 

.77980 

.75394 

.02586 

1.15226E-3 

30 

260 

500 

.69492 

.66327 

.03165 

1.43556E-3 

15 

300 

10 

.70192 

.68788 

.01404 

6.32160E-4 

15 

260 

50 

.59330 

.57749 

.01581 

5.87976E-4 

- 

0 

300 

.1 

.73014 

.72594 

.00420 

1.76310E-4 

0 

260 

1 

.53847 

.52832 

.01015 

2.40155E-4 

2500- 

45 

300 

1000000 

.84432 

.79373 

.05059 

4.56036E-3 

3000 

30 

300 

50000 

.75861 

.71130 

.04731 

2.7925  lE-3 

30 

260 

10000 

.87683 

.85817 

.01866 

5.74034E-4 

15 

300 

1000 

.72777 

.70690 

.02087 

8.53069E-4 
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300.01 


NEW  RADIANCE  ALGORITHM: 
APPROXIMATIONS  FOR  OPTICALLY  THICK  LAYERS 


Looking  only  at  a  single  isolated  layer,  straightforward  application  of  the 
simple  radiance  equation  leads  to: 

R  =  B  dr  or  R  =  (1-7)  B  (5) 

where:  R  =  Radiance, 

B  =  Planck  Function, 

and  dr  =  (1-7),  the  change  in  transmittance,  T,  across  the  layer. 

The  Planck  Function  is  defined  for  a  Curtis-Godson  density-weighted 
temperature  for  the  layer.  For  an  optically  thin  case,  the  observed  radiance 
in  this  simple  scenario  is  independent  of  viewing  direction. 

[Note  that  for  typical  lines-of-sight  across  a  multilayered  atmospheric  path, 
radiance  is  dependent  on  viewing  direction  while  total  transmittance  remains 
independent  of  the  observer’s  position.] 

However,  if  this  single  layer  is  optically  thick  and  includes  a  directional 
temperature  gradient,  the  observed  radiance  will  be  either  larger  (emanating 
from  a  warmer  thermal  region  closer  to  the  observer)  or  smaller  (emanating 
from  the  closer  cooler  region). 

A  simple  approximation  assumes  that  the  Planck  function  and  optical  depth 
vary  linearly  between  the  boundaries  of  the  layer: 

R  =  (1-7)  {  B„  +  2(B-B„)  [Ht  -  r/(l-7)]  }  (6) 

where:  n  implies  the  nearest  boundary, 

t  =  layer  optical  depth  between  boundaries, 
and  T  =  layer  transmittance,  again  for  a  single  layer. 
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NEW  RADIANCE  ALGORITHM:  Pg2 


For  a  multi-layer  scenario,  this  local  layer  radiance  becomes  the  input  for 
its  neighbor,  and  the  full  path  solution  is  reached  through  recursive 
calculations  and  the  layer  transmittance  then  becomes: 

r=  T(b  +  1)/T(b) 

where:  T(b)  and  T(b+ 1)  are  the  full  path  transmittances  from  observer 

to  boundary  b  and  b  +  1,  respectively. 

Because  of  the  degraded  2  cm'^  resolution,  the  required  optical  depth  term 
is,  in  reality,  an  "effective"  optical  depth,  also  derived  from  the  ratios  of 
adjacent  full  path  transmittances  at  MODTRAN  resolution: 

r  = -ln{T(b+l)/T(b)}.  (7) 

The  radiance  equation  for  a  single  layer  (Eq.  6),  including  the  "linear  in 
tau"  approximation,  appears  the  same  as  the  simple  radiance  equation  (Eq. 
5)  if  the  bracketed  quantity  {}  is  thought  to  contain  an  "equivalent"  Planck 
function,  defined  by  the  optical  depth  weighting. 
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COMPUTATIONAL  APPROACH 
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ANGLES  X  2  DIRECTIONS). 


MODTRAN  CALCULATIONS 


(qtu)  3ynss3bd 


COOLING  RATE  (K/day) 


MODTRAN  CALCULATIONS 


% 
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COOLING  RATE  (K/ciay) 


FASCODE  CALCULATIONS 
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Clough  et.  al.,  JGR,  97,  15761(1992). 


MAJOR  COMPUTATIONAL  DIFFERENCES 
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SUMMARY  OF  KEY  RESULTS 
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10.50  10.98  11.-17  11.95  12.44 

Wavelength(micron3)  Det  5  Filt  7(8  -  13  microns)(  769  -  1250  wn) 


SIMULATED  MEASUREMENT  CAPABILITY 
OF  THE  SPIRIT  III  INTERFEROMETER 

TANGENT  ALTITUDE  15  km;  RESOLUTION  2  WAVE  NUMBERS;  CHANNEL  10^13.0  MICRONS 
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SENSOR  NOISE 


SAMME  will  provide  a  code  ready  for  extensive  testing  in  the  CRAY 
environment. 

The  analysis  activities  will  include: 

1.  vector  optimization 

2.  instantaneous  full  spectral  analysis  of  heating/cooling 
rates 

a.  incorporation  of  molecular /thermal  climatological 
variability  (using  NRL-derived  climatologies)  and  assessment  of 
sensitivities  and  hemispheric  asymmetries 

b.  (as  above)  with  multiple  scattering 

c.  (as  above)  with  simple  cloud  types  and  new  NRL/DOE- 
derived  climatologies 

3.  calculations  of  photo-dissociation  rate  coefficients  with: 

a.  molecular /thermal  climatologies,  diurnal,  and  solar 
variability,  again  with  a  full  accounting  of  sensitivities  and 

hemispheric  asymmetries. 

b.  include  polar  day /night/twilight  scenarios  in  conjunction 
with  NRL/NASA. 

c.  (as  above)  with  assessment  of  uv-A/B  at  the  surface 

4.  exploration  of  energy  exchange  between  upper  and  lower 
atmospheres  as  a  function  of  climatological  and  auroral 
variability,  etc.  (NLTE) 

5.  implement  "alternate "  inversion  algorithm  schemes  for  any 
SEKDP,  ARM  or  other  fielded  instruments  with  compatible 
SAMME  spectral  resolution;  employ  compatible  simultaneous 
retrieval  algorithms. 
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SAMME  VIRTUES 


o  accomplish  with  one  radiance/irradiance  code: 

-  complete  with  solar  source  function, 

-  molecular  absorption-emission  descriptions, 

-  cloud,  aerosol,  and  molecular  default  profiles, 

-  validated  NLTE  rapid  algorithm 

-  plus  the  ability  to  input  new  data  as  appropriate 

o  determination  of  the  suitability  of  current  coarse  approximations: 

-  heating/cooling  rate  algorithms 

-  photodissociation  rate  coefficients 

-  energy  exchange  between  lower  and  upper  atmosphere 

o  current  climate-related  state-of-the-  art: 

-  extensive  spectral  degradation  for  efficiency 

-  UV  CFC  photodissociation  can  be  in  error  by  10-15% 

by  ignoring  Schumann-Runge  02  bands  T-dependence 

-  small,  systematic  errors  a  function  of  season  and  latitude 

-  similar  approximations  in  the  IR 

-  no  systematic  NLTE  coupling  with  lower  atmosphere 

o  exploration  of  this  net  impact  can  be  addressed  with  SAMME 

-  preliminary  assessement  this  summer:  IR  WMO/ICRCCM 
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CONCLUSIONS 


The  capabilities  of  the  MODTRAN2  code  have  yet  to  be  fully  exploited. 

o  Increased  accuracy  of  the  new  band  model,  coupled  with: 

-  the  LOWTRAN  7  and  FASCOD3  common  elements 

coase  continua  (CO2,  H2O,  N2,  O2,  etc.), 
spherical  refractive  geometry, 
default  constituent  profiles  for  gases,  clouds, 
aerosols,  fogs,  rains,  etc., 
molecular  and  particular  multiple  scattering, 
plus  ease  of  use 

o  MODTRAN2  may  be  effectively  employed  for: 

-  atmospheric  remote  sensing,  radiative  transfer 

0  MODTRAN2  transmittance  calculations  are  withing  %’s  of  FASCOD3 

-  statistically  and  in  spectral  detail 

-  for  simulations  at  2  cm  *  and  greater  resolution, 

MODTRAN2  may  be  substituted  for  FASCOD3 

-  layer-specific  radiance  contributions  are  excellent,  e.g.: 

detailed  agreement  in  the  Jacobian  comparisons 
primitive  cooling  rate  calculations  excellent 

o  MODTRAN2  flux-divergence  quantities  will  be  further  explored: 

-  optimized  estimates  of  up-  and  down-welling  fluxes 

-  appropriate  heating/cooling  rates  and  photodissociation  rates 

-  more  comparisons  with  experimental  data 

However,  if,  as  expected,  both  the  "total path"  and  "layer" 
integrity  of  the  MODTRAN2  calculations  are  maintained,  the 
community-at-large  should  find  this  new  tool  a  welcome  addition. 
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Retrieval  of  Tropospheric  Profiles  from  ER  Emission  Spectra: 
Preliminary  Results  with  the  DBIS 


J.-M.  Th^iault 

DREV-Defence  Research  Establishment  ValcartiCT 
P.O.  Box  8800,  Quebec,  Canada,  GOA  IRO 

G.  P.  Anderson,  J.  H.  Chetwynd,  E.  Mi^hy,  V.  Turner, 

M.  Cloutier  and  A.  Smith 

Geophysics  Directorate  /  PL,  Hanscom  AFB,  MA  01731  USA 
J.-L  Moncet 

Atmospheric  and  Environmental  Research,  Inc.,  Cambridge,  MA  02139  USA 


1.  INTRODUCTION 

Recently,  Smith  and  collaborators^’  ^  from  University  of  Wisconsin-  Madison  have  clearly  established  the  possibilities 
of  sounding  tropospheric  temperature  and  water  vapor  profiles  with  a  ground-based  uplooking  interferometer.  With  the 
same  perspective  but  for  somewhat  different  applications,  the  Defence  Research  Establishment  Valcartier  pi^V)  has 
initiated  a  project  with  the  aim  of  exploring  the  many  possible  avenues  of  similar  approaches.  The  central  objective  is  the 
development  of  methods  for  the  remote  sensing  of  atmospheric  profiles,  mainly  temperature  and  water  vapor,  that  affect  IR 
propagation  and  degrade  Electo-Optical  (EO)  system  performances.  There  are  several  important  issues  that  remain  to  be 
addressed  prior  to  the  definitive  implementation  of  such  a  technique  into  an  operational  system  that  responds  to  our  needs. 
In  order  to  address  some  of  these  issues,  DREV  in  collaboration  with  BOMEM  (Quebec  Canada),  has  developed  an 
instrument  referred  to  as  the  Double  Beam  Interferometer  Sounder  (DBIS).  This  sounder  has  been  conceived  to  match  the 
needs  encountered  in  many  remote  sensing  scenarios:  slant  path  capability,  small  field  of  view,  very  wide  spectral  coverage 
and  high  spectral  resolution.  Preliminary  tests  with  the  DBIS  have  shown  sufficient  accuracy  for  remote  sensing 
applications^.  In  a  series  of  field  measurements,  jointly  organized  by  the  Geophysics  Directorate  /  PL,  Hanscom  AFB,  and 
DREV,  the  instrument  has  been  run  in  a  wide  variety  of  sky  conditions.  Several  atmospheric  emission  spectra  recorded 
with  the  sounder  have  been  compared  to  calculations^’  ^  with  FASCODE  and  MODTRAN  models.  The  quality  of 
measurement-  model  comparisons  has  prompted  the  development  of  an  inversion  algorithm  ba^  on  these  codes.  The 
purpose  of  this  paper  is  to  report  the  recent  progress  achieved  in  this  research.  First,  the  design  and  operation  of  the 
instrument  are  reviewed.  Second,  recent  field  measurements  of  atmospheric  emission  spectra  are  analyzed  and  compared  m 
models  predictions.  Finally,  the  simultaneous  retrieval  approach  selected  for  the  inversion  of  DBIS  spectra  to  obtain 
temperature  and  water  vap)or  profiles  is  described  and  preliminary  results  are  presented. 

2.  INSTRUMENT  DESCRIPTION 

Essentially  the  DBIS  is  made  of  one  or  optionally  two  10-in.  diameter  Cassegrain  telescopes,  optically  coupled  to  a 
double-input  port  Fourier  transform  spectrometer  and  two  detection  units  (output  optics  1  and  2),  Figure  1  summarizes  the 
design  of  the  instrument.  This  configuration  allows  measurements  of  calibrated  spectra  according  to  the  following 
specifications:  any  selectable  zenith  angle,  scene  field  of  view  of  5  mrad,  spectral  coverage  from  3  to  20  |jm  and  a  spectral 
resolution  of  1  cm-1  or  greater.  As  a  part  of  the  input  modules,  a  large  flat  plate  scene  mirror  placed  in  front  of  each 
telescope  can  be  either  rotated  to  the  selected  scene  or  oriented  in  a  position  for  the  acquisition  of  calibrated  reference  spectra. 
The  pointing  capability  of  this  scene  mirror  allows  slant  path  measurements  from  0  to  360  degrees  with  a  tilt  adjustment  of 
±  10  degrees  in  azimuth  and  an  accuracy  of  0.1  degree.  The  coarse  adjustment  in  azimuth  is  simply  achieved  by  rotating  the 
whole  assembly  mounted  on  a  tripod.  After  reflection  on  the  scene  mirror,  the  beam  is  then  successively  focussed  by  the 
Cassegrain  telescope  and  reflected  by  an  off-axis  parabolic  mirror  (PM)  to  produce  a  collimated  b^  of  proper  diameter  at 
the  entrance  of  the  spectrometer.  The  two  output  modules  are  identical,  except  that  channel- 1  includes  a  MCT  detector 
optimized  for  the  5-20  |im  spectral  region  while  channel-2  includes  an  InSb  detector  optimized  for  the  2-5  pm  region.  Note 
that  the  InSb  module  is  not  explicitly  shown  in  Fig.  1.  These  modules  contain  parabolic  and  condensing  mirrors  that  focus 
the  beam  coming  from  the  interferometer  onto  a  detector  (MCT  or  InSb)  of  1  mm  diameter.  An  aperture  wheel  (AW) 
mounted  with  stops  of  different  diameters  permits  the  adjustment  of  the  field  of  view  of  the  instrument  (5  mrad  or  smaller). 
A  CCD  camera  integrated  into  the  telescope  module  can  be  used  to  aim  and  visualize  the  scene  under  consideration. 


192 


DOUBLE  BEAM  INTERFEROMETER  SOUNDER  (DBIS) 
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Figure  1:  DBIS  Optical  Train.  Figure  2;  Vertical  Profiles  of  T  and  RH 


Figure  3.  Atmospheric  uplooking  emission  spectra  as  measured  by  the  DBIS  system  and  calculated  by  FASCOD3. 


Figure  4.  Atmospheric  uplooking  emission  spectra  as  measured  by  the  DBIS  system  and  calculated  by  MODTRAN2. 
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INVERSION  ALGEBRA: 


Theriault  and  Moncet  (TM)  have  established  the  development  of  a 
successful  simultaneous  (temperature  and  water  vapor  profile)  retrieval 
algorithm,  based  primarily  on  FASCOD3  forward  calculations,  with 
accompanying  derivative  matrices. 

Traditionally  the  derivative  matrices  required  for  the  least  square  residual 
technique  embody  time-consuming  forward  runs  of  full-path  FASCODE 
radiance  predictions,  each  run  differing  from  its  predecessor  by  a  single  small 
perturbation, 

X  =  Xq  -f  x’, 

where:  x  =  T(K)  or  H20(g/m^),  for  each  layer,  /. 

The  Jacobian  matrix  is  then  defined  as  the  set  of  differences  in  total  radiance: 

dR(x,/)  =  R(x,Ib..Ro  (8) 

dx  x’ 

where:  R^,  is  the  unperturbed  total  radiance 

and  R(x,/)  is  the  total  radiance  with  a  single  perturbation 
(x  =  Xq  -t-  x’  and  x’  =  T’  or  H2O’)  at  layer  / 

The  size  of  the  original  matrix  is:  j  x  k 

where:  j  is  the  number  of  spectral  channels, 

often  dependent  on  spectral  resolution, 
and  k  is  (at  minimum)  the  number  of  atmospheric  layers 
or  boundaries  times  the  number  of  constituents 
undergoing  perturbation  in  the  simultaneous  retrieval. 
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M0DTRAN2  vs.  FASC0D3  JACOBIANS 


Moncet  and  colleagues^^  have  recently  devised  a  scheme  to  greatly  optimize 
calculations  of  the  Jacobian  elements,  based  on  FASCODE. 

However,  even  with  these  modifications,  the  task  still  consumes  a 
formidable  amount  of  computer  time. 

MODTRAN2  Jacobian  calculations:  each  full  path  radiance  calculation  has 
been  done  with  and  without  the  perturbation  at  each  layer  over  the  spectral 
range  of  the  DBIS  instrument. 

The  subsequent  derivative  matrix  elements  have  then  been  compared  to  the 
equivalent  FASCODE  elements. 

The  RMS  differences  in  the  Jacobian  radiances  (Eq.  8  with  the  denominator 
set  to  unity)  are  of  the  order  of  l.e-8  to  l.e-10  compared  to  an  average 
radiance  of  3.e-6  W/(cm2-ster-cm-l),  smaller  than  3  parts  in  1000. 
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Figure  5.  Temperature  Jacobians  for  a  2  K  temperature  perturbation  at  an  altitude  of  0.2  km  and  where  the  original 
temperature  and  water  vapor  profiles  correspond  to  supporting  radiosonde  data  appearing  in  Fig.  2  (jl02). 


Figure  6.  Water  vapor  Jacobians  for  a  O.lg  /  m^  perturbation  at  0.2  km  for  the  same  conditions  as  Fig.  5. 
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0.1  1  10 
ALTITUDE  (km) 


Figure  7.  Simultaneous  retrieval  of  T  and  H2O  profiles  for  the  hot-dry  case :  Inversion  from  SIMULATED  spectrum. 


ALTITUDE  (km) 

Figure  8.  Simultaneous  retrieval  of  T  and  H2O  profiles  for  the  hot-dry  case :  Inversion  from  MEASURED  spectrum. 
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PCModTRAN  2: 


ONTAR’S  PC  COMPATIBLE  MODTRAN  2  SOFTWARE 


L.W.  Abreu,  J.  Schroeder,  A.  McCann, 
J.  Kristi,  S.  Harvey  and  M.  Voltaire 


ONTAR  Corp. 

129  University  Road 
Brookline,  MA  02146 


PCModTRAN  2  is  an  implementation  of  the  Phillips  Laboratory/Geophysics  Directorate’s 
MODTRAN  2  model  and  contains  user-friendly  software  for  manipulation  of  the  input  and 
output  of  the  calculations.  This  package  is  compatible  with  IBM  and  all  other  personal 
computers.  The  software  package  has  an  input  generation  shell,  on-line  help  capability,  an 
ASCII  text  file  viewer  for  all  output  files,  screen  graphics  and  hard  copy  graphics.  A 
Cooperative  Research  and  Development  Agreement  for  MODTRAN  2  is  awaiting  final  approval 
by  Phillips  Laboratory. 
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PCModTRAN  2: 
ONTAR's  PC  COMPATIBLE 
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PCTRAN  7 


Personal  Computer  Version  of  the 
LOWTRAN  7  Atmospheric  Model 

Version  2 

June  1990 


ONTAR  Corporation 
129  University  Road 
Brookline,  MA  02146 
(617)-739-6607 


Air  Force  Geophysics  Laboratory 
Hanscom  AFB,  MA  01731 
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Run  n  1  of  1  MODTRAN  1  ^Nex-gr  <P1o^  <Exi^  <^10  =  Key^ 


Aeroaol  Model  Used  Navy  Maritime 


c 

T* 

o: 

r' 

r-  a 

(T  CQ 


X 

o 

w 

c 

-p 

.H 

2: 

cd 

u 

cr; 

> 

-p 

r 

0 

X 

X 

o 

o 

o 

o 

o 

cc 

CQ 

o 

o 

o 

o 

o 

c 

s 

X 

o 

o 

o 

o 

o 

-r" 

0 

P 

• 

* 

• 

■ 

9 

0 

t2l] 

c 

rH 

(T 

o 

4^: 

13 

X 

cd 

o 

c 

c 

:r3 

:z 

cc 

c 


X 

< 

0 

a 

y— V. 

-H 

a 

r" 

X 

a 

s' 

r- 

w 

S 

Ci] 

-H 

© 

© 

rH 

Q 

P‘ 

a 

a 

s 

w 

0 

o 

-H 

c 

p 

H 

-H 

D] 

T— {  ■ 

0 

w 

© 

0 

H 

H 

0 

1 

© 

-H 

> 

a 

© 

o 

a 

rH 

© 

o  ■ 

P 

> 

03 

CO 

c 

O 

P 

Cd 

© 

< 

© 

a 

© 

S 

P 

> 

p 

0 

>» 

< 

0 

a 

cd 

X 

p 

p 

Cd 

pi 

rH 

x 

w 

© 

d 

© 

> 

© 

0 

< 

X 

cd 

CO 

a  . 

a 

r. 

iH 

c 

H 

2: 

i«i^ 

0 

0 

0 

p 

p 

H 

1 

© 

0 

Pi 

a 

0 

o 

H 

X 

> 

•H 

© 

0 

p 

X 

P 

© 

o 

pi 

< 

p 

Cd 

© 

X 

£T3 

© 

© 

p 

S 

ft 

© 

0 

© 

p) 

< 

o 

01 

iM 

•H 

o 

00 

p 

© 

© 

cd 

c 

> 

> 

X 

\ 

X 

•H 

p: 

•H 

W' 

© 

cd 

c: 

s 

X 

ft 

cd 

© 

cm 

z 

•fp 

5 

pi 

0 

© 

K 

>5 

c 

! 

w 

-tH 

s 

0 

o 

\ 

cd 

X 

H 

E 

X 

p 

Pd 

© 

© 

© 

iH 

pH 

pi 

a 

p 

< 

© 

> 

< 

fd 

< 

a 

0 

© 

ft 

< 

Cd 

p 

cd 

rP 

p 

o 

CO 

K 

X 

0 

U 

s 

o 

o 

cd 

p 

c 

m 

© 

p 

X 

r* 

13 

nJ 

ft 

IP 

© 

© 

p 

c 

t 

H 

0 

© 

ft 

-H 

a 

a 

H 

Cd 

p 

CO 

XI 

< 

X 

X 

CO 

3r 

CN3 

Pd 

o 

209 


un  It  1  of  1  MODTRAN  Card  2  <Prev>  <Next>  <Hoine>  <Exit>  <F10  =  Keys> 


Initial  Altitude  (km)  _ qqO 

Final  Altitude/Tangent  Height  (km)  .000 
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Run  1  of  1  ONPLT  Scaling  Plot  11  1  FI f)=:Key^ 
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Run  #  2  of  4  LOWTRAN7  Cards  3  &  4 
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SAMM:  SHARC  AND  MODTRAN  MERGED 


A.  Berk,‘  D.C.  Robertson,*  L.S.  Bernstein,* 
R.L.  Sundberg,*  R.D.  Sharma,^  G.P.  Anderson,^ 
J.H.  Chetwynd,^  M.L.  Hoke,^  and  RJ.  Healey^ 


*Spectral  Sciences,  Inc.  Phillips  Laboratory/ GPOS 

99  South  Bedford  Street,  #7  29  Randolph  Road 

Burlington,  MA  01803-5169  Hanscont  AFB,  MA  01731-3010 

^Yap  Analytics 
594  Merrett  Road 
Lexington,  MA  02173 


Radiation  transport  calculations  of  ambient  IR  radiation  can  be  roughly  divided  into  two 
categories,  the  lower  altitude  regions  below  about  30  km  which  are  in  local  thermodynamic 
equilibrium  and  the  non-LTE  contributions  at  higher  altitudes.  MODTRAN  is  a  low  ^titude 
radiance  code  which  models  aerosol  and  molecular  scattering  and  absorptions,  refiractive  path 
geometry,  and  multiple  scattering.  SHARC  models  high  altitude  line-of-sight  radiances  by 
solving  the  pertinent  chemical  kinetic  equations  and  performing  equivalent  width  line-by-line 
(LBL)  calculations.  We  discuss  the  development  of  a  merged  model  which  has  been  dubbed 
SAMM,  SHARC  And  MODTRAN  Merged.  SAMM  includes  a  fast  expanded  LBL  approach 
which  is  applicable  to  both  altitude  regimes,  and  contains  all  MODTRAN  and  SHARC  model 
features.  We  will  demonstrate  that  SAMM  predictions  vary  smoothly  as  LOS  tangent  heights 
are  increased  from  low  to  high  altitude. 


214 


SAMM:  SHARC  AND  MODTRAN  MERGED 
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MOTIVATION 
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SAMM  IS  THE  SOLUTION,  A  SINGLE 
CORRELATED  LOS  RADIANCE  MODEL 
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CALCULATIONAL  SEQUENCE 
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INPUT  MODULE  (SAMPLE  MENU) 
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NON-LTE  POPULATIONS  MODULE 
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LAYER  -  TO  -  LAYER  PUMPING 


REFRACTIVE  GEOMETRY  MODULE 
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DETERMINES  OPTICAL  AND  SOLAR 
PATH  COLUMN  DENSITIES 


RADIATION  TRANSPORT  MODULE 
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2-FLUX  MULTIPLE  SCATTERING  CALCULATION 


INITIAL  VALIDATION 
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SOME  ILLUSTRATIVE  VALIDATION  CASES 
ARE  PRESENTED. 


MIDDAY 

LIMB,  TANGENT  TO  GROUND 
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WAVELENGTH  (jim) 
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SCATTERING  CAN  DOMINATE  HIGH 
ALTITUDE  LOS  RADIANCES  BELOW  4 //m. 


THE  MODERATE  SPECTRAL  ATMOSPHERIC  RADIANCE 
AND  TRANSMITTANCE  (MOSART)  PROGRAM 


W.M.  Coraette  D.C.  Robertson 


Photon  Research  Assoc.,  Inc. 
10350  N.  Torrey  Pines  Road 
Suite  300 

La  Jolla,  CA  92037-1020 


Spectral  Sciences,  Inc. 

99  South  Bedford  Street,  #7 
Burlington,  MA  01803-5169 


The  MOSART  program,  a  new  computer  program  for  predicting  and  evaluating  the  radiative 
environment,  incorporates  features  from  the  MODTRAN  2  and  APART  (Version  7.00) 
programs,  plus  adding  some  additional  new  features.  The  MODTRAN  2  1  cm'’  and 
LOWTRAN  7  20  CM'*  band  parameters  and  geometry  specifications  are  combined  with  the 
APART  global  climatology  data  bases,  global  terrain  background  data  bases,  and  ray  tracing 
algorithm  to  produce  a  combined  atmosphere  and  background  program  for  supporting  target  and 
background  scene  modelling,  in  particular  the  Strategic  Scene  Generation  Model.  The  program 
can  interpolate  primary  atmospheric  parameters  over  multiple  model  atmospheres  as  well  as  use 
a  single  model  atmosphere.  Heat  transfer  calculations  for  the  terrain  materials  are  included,  and 
terrain  clutter  parameters  (i.e. ,  mean,  standard  deviation,  PSD)  are  output.  Utilities  for  creating 
input  files,  spectral  tables,  spectral  plots,  luminance,  and  statistical  two-dimensional  backgrounds 
are  included. 
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Dr.  David  C.  Robertson 
Spectral  Sciences,  Inc. 
Burlington,  Massachusetts 
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strategic  Scene  Generation  Modei  Requirements 


SSGM  Requirements  lor  MOSAR  F 
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Broad-Band  Thermal  and  Solar  Loading 
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THERMAL  PATH  RADIANCE 
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SCATTERING  ANGI  E  (dcol 


1  Jne  Correlation  Eflects 
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Photon  Research  Associates.  Inc. 


WAVELENGTH  (  iim) 
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Photon  Research  Associates,  Inc. 


LU 

C/) 


CO 


CM 

C 

CQ 

o 


04 

o 


O 

c  ■*- 

O  CM 


(0 

C/) 

■  ■■■ 

E 

111 

15 

E 

'imm 

o 

1- 


O 

r 

o 

> 

c 

o 

O 


U) 

c 

•w 

im. 

0) 

fS 

o 

(/) 

O) 

'o 

ro 

CC 


o 

CO 

D 

+ 

i 

CD 

CD 

CM 

1 

+ 

CM 

c 

CM 

c 

03 

CM 


II 

CO 


o 

IT 

(O 

04 

o 


04 

(0 

o 


(A 

O 

CO 

3 

O 

3 

LL 

>> 

'w 


0) 

CO 

IT 

(0 

rj 

04 

o 


S'  o 


o 

Q 

o 

.n 

E 

3 


arcM 

.—  Q 

C  2 

5  o 

Q  2 

CM 

<  e- 

^  o 


CM 

■  ■ 

o 

c 

CM 

C 

V) 

o 

C 

O 

c 

o 

im. 

■■w 

CO 

O 

r 

o 

o 

o 

< 

3 

■c 

0) 

> 

c 

o 

#ii< 

CO 

o 

CO 

o 

E 

o 

3 

LL 

> 

c 

o 

o 

3 

Vi 

£ 

o 

o 

o 

V) 

•  ■M 

o 

CO 

(/) 

o 

< 

CO 

1 

o 

1 

1 

249 


Atiiiospheric  Characterization 
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ATMOSPHERIC  CHARACTERIZATION 
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Cirrus  Clouds: 

-  Standard  (64  nm  +  Extinction) 

-  Subvisual  (4  jim  +  Extinction) 

-  Heymsfield  (Temperature  Dependent) 
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Ion  Hetearcn  Astociaies,  me.  _ _  m.  m  m.  m  « ■  m  m  « 

GLOBAL  CLIMATOLOGY  PROFILE:  11  jim;  NADIR  VIEW; 

0°  LONGITUDE;  21  JUNE  1993,  12:00 
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Mean  Star  Radiance 
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TROPICAL  ATMOSPHERE 


AIR  TEMPERATURE 
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SCNGEN:  Create  Deterministic/Stalistical  Terrain  Scenes 
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B  OlS  tl.lS 


INCLUSION  OF  ACCURATE  MULTIPLE  SCATTERING  IN  MODTRAN 


K.  Stamnes  and  S.  Tsay 


Geophysical  Institute 
Univ.  of  Alaska 
Fairbanks,  AK  99775-0800 


M.  Yeh 

Caelum  Research  Corp. 

Silver  Spring,  MD  20901-4554 


Our  work  on  the  MODTRAN  computational  code  aimed  at  improving  the  computation  of 
multiple  scattering  not  only  in  the  infrared  part  of  the  spectrum,  but  also  in  die  solar  part  where 
scattering  effects  are  important  under  both  cloudy  and  clear  sky  conditions.  This  task  is 
concerned  with  the  inclusion  of  the  general-purpose  multiple  scattering  algorithm  DISTORT 
(Discrete  Ordinate  Radiative  Transfer  code  as  summarized  by  Stamnes  et  al.,  1988)  into 
MODTRAN,  including  proper  interfaces  to  allow  for  accurate  computation  of  multiple  scattering 
effects.  In  the  solar  part  of  the  spectrum  this  is  important  for  clear  sky  conditions  (Rayleigh 
scatter)  as  well  as  cloudy  and  hazy  situations.  In  the  terrestrial  infrared  part  of  the  spectrum 
(beyond  4  microns)  molecular  (Raleigh)  scattering  is  negligible,  but  scattering  from  clouds  and 
aerosols  can  not  be  ignored.  Results  were  compared  between  the  two  versions  and  their 
differences  and  improvements  will  be  discussed. 
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Inclusion  of  Hccurate  Multiple  Scattering  in 

MODTRHN 
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»  Multiple  scattering  in  MODTRBN  is  based  on  a  2-streani  code  (BMFLUH) 
uiith  an  isothermal  layer  appronimation  from  mhich  iipmard  and  domnuiard 
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accuracy  of  retrieuals  of  remotely-sensed  atmospheric  properties  that  rely 
the  interpretation  of  measured  radiances  (yround-based  or  from  space). 


do  we  Improve  Multiple  Scattering  (MS) 

Treatment? 
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*  extensive  testing  is  required  und  must  be  carefully  executed.  This  is  a  time- 
consuming  undertaking. 


Figure  1.  Rowchart. 
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Radiance  (W.m'^.sf'41  m'^)  Radiance  (W.m’  .sr'  4im‘ )  Radiance  (W.m"  .sr 


Wavelength  ( p,  m) 


Fi«nire  2.  Comparison  of  solar  multiple  scattering  source  function  for  a  dear  sky  su^ 
°  arctic  winter  McClatchey  atmosphere,  layer  1  for  the  lowest  atmospheric 
layer  and  layer  32  for  the  highest  layer. 
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Wavelength  (  um) 


Figure  4.  Same  as  in  Figure  3,  but  for  vertical  profiles  of  source  functions  at  0.6  |im 
and  spectral  total  radiance  at  top  for  an  observer  looking  straight  down. 
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Atmospheric  Layering 
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uery  efficient.  Therefore  it  uiould  be  of  great  interest  to 
find  out  under  uihat  conditions  it  is  uaiid. 


B.  Multiple  Scattering  in  Sphericai  Geometry 
Possible  approach: 

(i)  Use  plane  geometry  to  approKimate  deriuatiue  term,  but  compute  Chapman  function 
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2.  For  Isotropic  scattering  there  is  nothing  ‘driuing*  azimuth-dependence.  So  solution  mill  be  azimuth- 
independent  in  this  case. 


.  Multiple  Uieiuing  Directions 


E 

o 

Im 

W 

C/) 

1. 

(U 

ID 

C 

CD 


u 

0) 

L. 

9mm 

•a 

05  ^ 

£  S 

3  i: 

O  C 

•-  o 
3  u 

03  C 

05  . 

c  *a 
—  0) 

fO 

03 

L.  *0 

°  .2 

(/)  c 
0}  a 

3  u 
ca  03 

o  ^ 

“  3 

35  03 

;=  = 

c  re 

03  c 

s| 

a  3 


a 


re 

L. 

U 

(/) 


re 

u 


c 

05 

'3) 

C 


re 

50 

o 

.2 

.£ 


3 

s 

3 

3 

■3 

.£ 

50 

3 

■3 

3 

U 

3 

•3 

3 

3 

3 

? 

2 

s 

re 

3 

re 


u  c 


c 

3 


3 

CO 

O 

£ 

3 

J3 

2 

3 

O 

3 


3 

U 

*3! 

3 

■o 

3 

L. 

re 


CO 

3 


3 

O. 

O 


C/3 

G 

c 

u 

G 

IMS 

ISM 

Q. 

u 

O 

(/) 

G 

o 

G 

u 

‘"5 

3 

75 

ra 

a 

c 

o 
■  ■■■ 

#1^ 

G 

IM« 

£ 

£ 

o 

u 

re 

3 

•  MB 

o 

3 

3 

Q. 

3 

Z 

£ 

#11% 

oL 

iMi 

s 

CC 

Q 

4-^ 

H- 

u 

3 

75 

c 

£ 

o 

i: 

o 

•■■■ 

IMS 

3 

CO 

IMS 

■a 

C/) 

3 

3 

“O 

CO 

w 

CO 

u 

3 


3 

C 

O 


u 

re 


3 

1^ 


c 

re 

3 

u 


re  re 

3  3 
U  31 

_  c 
^  re 

1 1 

£  „ 

g,“ 

05  3 

re  o 

^  CO 

re  ? 

a.E 

35  u 

^  Si 
«  re 
£  “ 
o 

0)  3 
05  05 

CO  3 

Z 

C  3 
CC  £ 


C 

o 

re 

3 

O. 

£ 

o 

u 

3 

u 

C 

.2 

■o 

re 

L. 

3 

3 

a 

£ 

o 

u 

3 


O 

X3 


I— 

CC 

a 

C/l 


3 


O 

o  ^ 

"  I 

u 


o 


o 

CO 

U 

3 

3 

£ 

3 

C 

3 

3 


303 


Accurate  Computation  of  Heating/Cooling 
Rates  and  Potoigsis  Rates 
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RAYLEIGH  AND  AEROSOL  SCATTERING  IN  THE  TROPOSPHERE  AND 
STRATOSPHERE  IN  THE  SPECTRAL  RANGE  175-850nin:  AN  INTERACTIVE  MODEL 


D.E.  Anderson 


R.  De  Majistre  and  S.  Evans 


The  Johns  Hopkins  Univ.  Computational  Physics,  Inc. 

Applied  Physics  Lab  Fairfax,  V A  22231 

Laurel,  MD  20723 


We  have  developed  an  interactive  version  of  the  radiation  field  model  described  by  Anderson 
and  Lloyd  (1990),  Anderson  and  DeMajistre  (1992)  and  DeMajistre,  Zasadil  and  Anderson 
(1992).  The  model  provides  solution  for  the  diffuse  radiation  field  including  ground  albedo, 
aerosols  and  absorption  by  molecular  oxygen  and  ozone.  The  integral  equation  solution  has  been 
coded  for  representation  of  the  radiation  field  and  accurate  (—  10  ->  20%)  approximations  to 
anisotropic  scattering  are  included.  LOWTRAN  model  atmospheres  and  Rayleigh,  aerosol  and 
ozone  optical  parameters  have  been  utilized  in  the  model.  A  spherical  shell  model  for  initial 
energy  deposition  is  employed. 
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Altitude  (km)  Altitude  (km) 
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Aerosol  Extinction  (volcanic) 


Altitude  (km)  Altitude  (km) 
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Plot 


We  have  used  the  model  to  calculate 
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EXPERT  SYSTEM  IN  PLEXUS 


P.C.F.  Ip,  S.B.  Downer,  M.  Noah,  F.O.  Clark 

K.  Radermacher,  J.P.  Kennealy 


Mission  Research  Corp. 
One  Tara  Blvd. ,  Suite  302 
Nashua,  NH  03062 


Geophysics  Directorate, 

Phillips  Laboratory 
Hanscom  AFB,  MA  01731-3010 


This  talk  will  review  the  principles  behind  building  the  PLEXUS  expert  system.  The  use  of 
an  expert  system  in  the  PLEXUS  interface  will  intelligendy  assist  users  in  solving  atmospheric 
background  problems  by  directing  them  to  the  appropriate  atmospheric  codes(s)  and  only 
querying  them  for  parameters  important  to  their  application.  Building  the  PLEXUS  expert 
system  is  a  long  and  detailed  process  involving  several  iterative  phases.  The  knowledge  base 
is  acquired  by  researching  existing  code  documentation,  consulting  experts  in  the  field,  and 
establishing  trends  from  running  the  atmospheric  codes.  As  the  knowledge  base  is  developed, 
it  is  organized  into  meaningful  subtopics  and  then  translated  into  rules  used  by  the  expert  system 
shell.  Examples  of  each  of  these  phases  will  be  given,  with  main  emphasis  placed  on  the  results 
obtained  from  a  systematic  study  of  MODTRAN  input  parameters.  TTiis  study  has  yielded  some 
results  of  interest  to  a  large  number  of  MODTRAN  users. 
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KO.  Clark 

Geophysics  Directorate,  Phillips  Laboratory 

Hanscom  A  FB,  MA 


PLEXUS  Expert  System  Goals 
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PHILLIPS  LABORATORY’S  EXPERT-ASSISTED  USER  SOFTWARE 

(PLEXUS) 


* 


S.B.  Downer,  J.P.  Kennealy,  P.C.F.  Ip,  F.O.  Clark 
M.  Noah,  KJ.  Radermacher,  D.  Einstein 


Mission  Research  Corp. 
One  Tara  Blvd.,  Suite  302 
Nashua,  NH  03062 


Geophysics  Directorate, 

Phillips  Laboratory 
Hanscom  AFB,  MA  01731-3010 


PLEXUS  Version  1.0  is  now  available  for  general  release.  This  version  supports  the  DoD 
Atmospheric  codes  LOWTRAN,  MODTRAN,  SHARC,  the  CBSD  Celestial  Codes  and  the 
FAUST  Method.  The  FAUST  Method  is  designed  to  help  fulfill  current  DoD  requirements  for 
downward-looking  background  radiance  and  transmittance  descriptors.  By  combining  the  results 
from  the  equilibrium  lower  atmospheric  model  (MODTRAI^  and  non-equilibrium  upper 
atmospheric  model  (SHARC),  a  unified  seamless  single  answer  is  generated.  The  FAUST  user 
interface,  coupled  to  an  intelligent  knowledge-based,  assists  users  of  all  abilities  in  correctly 
setting  up  and  running  calculations  of  this  nature. 

PLEXUS  provides  a  number  of  data  visualization  tools  as  well  as  importing,  exporting,  and 
printing  functions.  Two  such  tools,  integrated  band  intensities  and  filter  convolutions  allow 
users  to  explore  sensor  operability  ranges.  Present  and  future  capabilities  are  presented  with  a 
demonstration  of  the  first  general  release,  PLEXUS  Version  1.0. 
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FAUST  Method  Interface 
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Correct  Code  Selection 
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EARTHLIMB  BACKGROUNDS  IN  THE  STRATEGIC  SCENE  GENERATOR  MODEL 


Susan  McKenzie 


Mission  Research  Corporation 
One  Tara  Blvd. ,  Suite  302 
Nashua,  NH  03062 


The  Strategic  Scene  Generator  Model  (SSGM)  is  a  computerized  methodology  which  generates 
viewer-perspective  radiance  maps  of  quiescent  and  enhanced  natural  backgrounds  and  perturbed 
backgrounds  with  embedded  targets  and  target  induced/related  events.  The  Earthlimb 
Backgrounds  represents  an  on-line  phenomenology  capability  for  the  SSGM  Baseline  design. 
The  components  for  this  phenomenology  have  been  provided  by  Phillips  Laboratory  codes. 
RAD_E  v4.4  is  the  computer  model  that  generates  structured  Earthlimb  2D  radiance  maps  in 
the  SSGM.  The  ID  radiance  profiles  are  computed  by  one  of  two  options:  on-line  calculations 
by  PL  radiance  codes  or  extractions  from  the  FASTLIMB  database.  FASTLIMB  is  a  1% 
resolution  database  generated  off-line  using  PL  radiance  codes  and  is  supplied  for  rapid-running 
profile  generation.  Both  options  store  in-band  ID  radiance  profiles  and  structure  information 
for  generating  2D  image  output.  2D  structured  images  can  be  generated  by  RAD_E  directly. 
This  structure  contains  stochastic,  non-homogeneous,  non-stationary  spatial  structure.  A 
correlated  multi-frame  option  is  included  in  this  process  for  frame-to-frame  continuity .  RAD_E 
radiance  and  structured  image  output  generated  using  these  methodologies  will  be  presented. 
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Two  Operational  Modes: 
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ATMOSPHERIC  EXTINCTION  ANALYSIS  FOR  THE  AIR  FORCE 
AIRBORNE  LASER  PROGRAM 


Larrene  K.  Harada  and  Daniel  H.  Leslie 


WJ.  Schafer  Associates 
Arlington,  VA 


The  Air  Force  is  developing  an  airborne  laser  system  capable  of  engaging  targets  at  ranges 
of  at  least  several  hundred  km,  at  elevation  angles  from  1  to  4  degrees.  We  will  present  our 
wavelength-trade  analysis  comparing  the  capability  of  several  candidate  laser  systems.  The 
impact  of  molecular  and  aerosol  scattering  and  absorption  has  been  assessed  using  FASCOD2 
and  the  HITRAN92  database,  and  these  results  will  be  presented.  We  are  using  the  NASA 
SAGE  satellite  data  base  to  provide  recent  stratospheric  extinction  data,  and  comparisons  with 
standard  LOWTRAN  models  will  be  described. 
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A  FAST  SCHEME  FOR  A  LINE-BY-LINE  FORWARD  MODEL 


S.  Kadokura,  A.  Shimota 


A  fast  scheme  for  a  line-by-line  forward  model  has  been  developed.  The  keypoint  of  the 
scheme  is  a  hybrid  of  analytical  and  numerical  integration,  and  it  reduces  the  calculation  time 
approximately  by  an  order  of  magnitude  as  compared  with  that  of  the  classical  scheme.  The 
scheme  is  summarized  as  follows:  (a)  divide  the  whole  range  of  the  wave  number,  v,  into  many 
intervals  with  various  widths,  so  that  a  quadric  is  obtained  as  an  approximation  of  the  cross 
section  for  each  interval  (making  a  data  base  of  the  cross  section);  the  transmittance  for  a  wave 
number  interval,  T,  is  the  integral  of  the  exponential  function  of  the  optical  depth,  ,  where 
,  is  expressed  approximately  with  a  quadric  by  that  approximation,  i.e.  (b)  if 

a>0,  calculate  the  exponential  or  error  function  to  obtain  the  integral;  (c)  if  a<0,  divide  the 
interval  into  sub-intervals  where  the  terms  of  v^  can  be  neglected,  so  that  the  transmittance  for 
the  sub-interval,  T  ,  is  expressed  with  the  exponential  function;  calculate  the  weighted  mean  of 
T  to  obtain  T. 
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APPLICATION  OF  FASCODE  PROGRAM  TO  A  HIGH  TEMPERATURE  GAS 
MONITORING 


N.  Takeuchi 


Remote  Sensing  &  Image  Res.  Center 
Chiba  University 

1-33,  Yayoi-cho,  Inage-ku,  Chiba-shi 
Chiba  263  Japan 


FASCODE  program  was  used  to  investigate  the  interference  absorption  by  CO2,  which  affects 
the  monitring  of  CO  concentration  at  a  high  temperature  state.  A  tunable  lead  salt  diode  laser 
was  used  in  4.5  micron  region  and  it  was  found  that  the  hot  band  gas  absortion  is  not  satisfactly 
explained  by  the  program. 
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Spectroscopy  of  Combustion. 

■  Flame 

■  Engine 

■  Exhaust  Flue 

Physical  Quantities  to  be  measured: 

■  Velocity 

•  Temperature 

•  Pressure 

■  Density 

•  Concentration  of  Chemical  Species 

•  Particles 

■  Size  Distribution  of  Fuel  Particles 
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Mix.  Ratio 


Current  Topics  of  Combustion. 


Intermediate  Product 


Flame  Temp. 

- . .  - . 

High 

Low 

Lower 

CO 

NOx 

Higher 

NOx 

CO 

CH4  ^  CHs  CH  +  N2  —  HCN 

— ►  •••  — ►  NOx 
I  400 °C  No  growth  of  reaction 
1000°C  C0  +  0H^C02  +  H 

preferable  spatially  homogeneous  temperature 
distribution  to  burn  at  higher  temperature. 
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WAVENUMBERS 


v’-v" 

Vo(cm-^) 

Vmin-Vmax 

IS 

Smax 

1-0 

2143.27 

1911-2288 

1.036E-17 

4.636E-19 

2-0 

4260.06 

3985-4361 

7.692E-20 

3.532E-21 

3-0 

6350.44 

6032-6418 

4.909E-22 

2.338E-23 

4-0 

8414.47 

8053-8465 

1.531E-24 

8.159E-26 

v'-v" 


10°02-01'01 

01’01-00°01 

02'01-01'01 

10°01-01'01 

02'11-02‘01 

10°11-10°01 

i0°n-10°02 

ol'n-m’m 

oo°n-oo°oi 

10°12-00°01 

10°11-00°01 


Vo(cm-i) 


618.03 

667.34 

667.75 


2324.14 
2326.60 
2327.43 
2336.63 

2349.14 
3612.84 
3714.78 


V  min-Vmax 

IS 

546-687 

1.46E-19 

593-752 

8.02E-18 

600-750 

6.09E-19 

649-791 

1.564E-19 

2227-2371 

2.601E-19 

2231-2372 

1.021E-19 

2231-2374 

1.716E-19 

2227-2384 

7.0nE-18 

2230-2397 

9.157E-17 

3509-3661 

1.003E-18 

3610-3763 

1.504E-18 

Smax 


5.459E-21 

2.982E-19 

1.122E-20 

5.850E-21 

4.985E-21 

3.929E-21 

6.621E-21 

1.344E-19 

3.526E-18 


3.855E-20 
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100,000 
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FASCODE  calculation  of  the  1-0  CO  absorption  band 


L<Itl 


CO,  CO2,  and  H2O  Absorption  Spectrum  Calculated 
with  FASCODE.  (  CO  SOppin,  CO2  15%,  H2O  20%  ) 
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Wavenumbers 

CO,  CO2,  and  H2O  Absorption  Spectrum  at  400 °C 

Calculated  with  FASCODE 
(  CO  SOppin,  CO2  15%,  H2O  20^ 


Absorption  at  2169.1984cm‘^  (Room  Temperature) 
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g(v)  dv 


im.t  un.t  tm.i  ttn,4  nai.t  titt  iim.i  iim.i  KM.t  iiN.i  im.i  tin.* 


MavanuMbtri  n,  c 


Room  Temperature  400°C 

Calculated  Spectrum  at  2169.2cm'l  with  FASCODE 


Room  Temperature  400°C 

Observed  Spectrum  at  2169.2cm’^ 
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experimental  error  ? 


